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Abstract. An ElectroMagnetic Solver (EMS) [Chen et al , Phys. Plasmas, 13, 
123507 (2006)] is employed to model a spectral gap and a gap eigenmode in a 
periodic structure in the whistler frequency range. A Radially Localised Helicon 
(RLH) mode [Breizman and Arefiev, Phys. Rev. Lett, 84, 3863 (2000)] is considered. 
We demonstrate that the computed gap frequency and gap width agree well with 
a theoretical analysis, and find a discrete eigenmode inside the gap by introducing 
a defect to the system's periodicity. The axial wavelength of the gap eigenmode is 
close to twice the system's periodicity, which is consistent with Bragg's law. Such an 
eigenmode could be excited by energetic electrons, similar to the excitation of Toroidal 
Alfven Eigenmodes (TAE) by energetic ions in tokamaks. 



Gap eigenmode of radially localised helicon waves 



2 



1. Introduction 

It is a generic phenomenon that spectral gaps are formed when waves propagate in 
periodic media, and eigenmodes can exist with frequencies inside the spectral gaps if a 
defect is introduced to break the system's perfect translational symmetry, and create 
an effective potential well to localise these waves. [H El El HI El [6] Fusion plasmas have 
a few periodicities that can produce spectral gaps: geodesic curvature of field lines, [7] 
elongation [HI |9] or triangularity of flux surfaces, [10] helicitypTj and periodic mirroring in 
stellerators.[12j Weakly damped eigenmodes which are readily destabilised by energetic 
ions often reside in these gaps, and they may degrade fast ion confinement. [T3l [T4] The 
most extensively studied gap eigenmode is the toroidcity-induced Alfven eigenmode 
(TAE),[15] however, there are also numerous other modes with similar features. [T6] 

Zhang et al [TT] observed a spectral gap in the shear Alfven wave continuum in 
experiments on LArge Plasma Device (LAPD) with a multiple magnetic mirror array, 
and obtained consistent results through two-dimensional numerical modelling using the 
finite difference code, ElectroMagnetic Solver (EMS).[18j Although eigenmodes inside 
this gap were not formed, a possible experimental implementation was proposed to 
detect them. The idea is to use the endplate of the machine as a "defect" that breaks 
the axial periodicity or to vary current in one of the independently powered magnetic 
coils. 

In the present study, we will use EMS to examine a gap eigenmode in a linear system 
with slightly broken axial periodicity. We will consider the Radially Localised Helicon 
(RLH) mode, whose radial structure has been described in [T9j . We will show that the 
computed gap frequency and gap width agree well with a theoretical analysis, and that 
there is a discrete eigenmode inside the gap. Such an eigenmode could be excited by 
energetic electrons, similarly to the excitation of Toroidal Alfven Eigenmodes (TAE) by 
energetic ions in tokamaks. 

2. Theoretical analysis 

2.1. Basic equations 

The spatial structure of RLH wave field and the corresponding dispersion relation can 
be found from (17) in [19] that applies to a whistler- type linear wave in a cold plasma 
cylinder and has the form: 

1 d dE m 2 g m uo 2 Edg/dr 

r dr dr r 2 k 2 r c 2 1 + (mdg / 'dr) jk 2 r^ 

where E = E z — (k z r /m)E Lp with E z and E^ the axial and azimuthal components of the 
wave electric field, respectively. Functions g and rj are g = uj 2 /ujuj c and rj = —oj 2 /oj 2 , 
respectively, where uj p is the electron plasma frequency and uj c the electron cyclotron 
frequency. 

In this section, we will limit our consideration to the case of sufficiently dense 
plasma in which uj p /c 3> m/a, where a is the plasma radius. This assumption makes it 
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allowable to set E z = and ignore the term (mdg/dr) /k^rrj in the denominator of (1), 
which gives 



We now generalise this equation to the case of slightly modulated (z-dependent) plasma 
equilibrium with a following separable form of (j0 2 /u) c ' 

— — [1 - e(z) cos gzj . (3) 

Here, 6«1 and g are the modulation envelope and wavenumber respectively. The axial 
scale-length of the envelope e(z) is assumed to be much greater than 1/q. We recall that 
k 2 z = —d 2 /dz 2 and transform (2) to 



h (I r | r ^- m ^) = ^[1 -e(z) cos q z]r-^ fa± 



(4) 



The modulated equilibrium introduces resonant coupling between the modes with 
k z = q/2 and k z = — g/2, which suggests the following form for 

E^ = A+e %qz l 2 + A_e- %qz l 2 , (5) 

where A + and A_ are slow functions of z compared to cos qz. 

Let ujq be an eigenfrequency of (4) for e = and |fc z | = q/2. It is then 
straightforward to separate spatial scales in (4) and obtain a set of coupled equations 
for A i and A_: 



-£(&r&rA_-mM_)-A_m$r|: (W 



(6) 



(7) 



We have neglected second axial derivatives of A+ and A_ and arranged (6) and 
(7) so that their right hand sides represent relatively small terms compared to the 
dominant terms on the left hand sides. With this ordering, we conclude that the radial 
dependencies of A+ and A_ need to be close to the eigenfunction \&(r) of the lowest 
order ODE 

.2 
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( |-r|-r* - m 2 *) + W=?r£ (^\ = 0. (8) 
\or or J c z or \cj c o J 

The differential operator in this equation is self- adjoint. As a result, multiplication of 
(6) and (7) by and integration over radius lead to 



= ~iqf z I drr* (&r&rA+ - m*A + ) + / drr* A + m^r ■& ( £ 

(9) 



-\ f drr^A_m<fe{z)r A 



pO 

Or \ cuco 



Gap eigenmode of radially localised helicon waves 

= iq§-J drr* (&r&r A. - ra^A.) + / drr^A.m^r^ 
-Udrr*A +m ^e(z)ri;(^). 



(10) 



Having eliminated the lowest order terms, we now set 

A + = F(z)V(r), A_ = G{z)^{r) (11) 
in (9) and (10) to obtain the following set of coupled equations for F(z) and G(z): 

= - — + F °- - -e(z), 12 

q oz ujq 2 

= — + G -^e^). (13) 

Spectral gap and continuum 
In the case of z-independent e, (12) and (13) have exponential solutions 

FocGoc e iKZ (14) 
with k the wave number. The corresponding dispersion relation has two roots: 



(15) 



= uo[l + V(4«/?) 2 + (e/2) 2 ], 

= w [l - ^(4«/g)2 + (€/2)2]. 

Here, clJ+(ac) and are the continuum frequencies above and below the spectral gap, 

respectively. The upper and lower tips of the spectral gap correspond to k = 0: 

"+(0)=o;o(l + §), 
cj_(0)=o;o(l-f). 

The normalised width of the gap, Auj = [u;+(0) — cj_(0)]/cl;o, is therefore equal to the 
modulation amplitude e. The central frequency of the gap ujq needs to be found from 
(8) as the RLH eigenfrequency for an axially uniform plasma cylinder and it can be 
written as 

c, c q(0)cV 

4^ (0) 2 ' V } 

where cj c0 (0) and u) p o(0) are the on- axis values and the numerical form-factor V is 
determined by the plasma radial profile. 



2.3. Wall-localised eigenmodes 

A discrete-spectrum eigenmode can be created inside the spectral gap if the system's 
periodicity is broken. We illustrate that by considering an ideally conducting endplate 
located at z = z , so that the plasma now occupies only a half cylinder to the right of 
the endplate [z > zq). The boundary condition at the endplate is 

E (p (r;z o ) = (18) 

or, equivalent ly, 

F(z )e iqz °/ 2 + G(z )e- iqz °/ 2 = (19) 
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(see (5) and (11)). The electric field of the discrete-spectrum mode must also vanish at 
z —> oo, i. e., 

F(oo) = G(oo) = 0. (20) 
Equations (12) and (13) (with e = const) admit an exponential solution 

FocGoc e~ Xz (21) 

with a positive value of A that satisfies boundary conditions in (19) and (20). 

Indeed, the exponential ansatz, (21), reduces (12), (13) and (19) to an algebraic set 



4iA _ u — uj{] G 

= F + F --ire, 

q uj 2 

q ll> 2 



(22) 
(23) 



F + Ge~ iqz ° = 0. (24) 
The solvability conditions for this set determine u and A as functions of zq as follows: 

(25) 



uj — ojo e 

--cos qz , 

uq 2 



A — -qsinqzQ. 



(26) 



We observe that A is positive for < qz < 7r, and the wave field therefore vanishes at 
infinity in this case to represent a discrete eigenmode whose frequency is automatically 
inside the gap. In particular, the eigenfrequency is exactly at the gap centre when 
qz = 7r/2. We further note that this wall-localised solution can also be viewed as an 
odd-parity eigenmode (E^(z) = — £^(2z — z)) in the entire periodic cylinder with a 
defect at z = zq. An axial profile of (jjI/(j0 c with such defect is shown in figure 1(a). 




(b) 




Figure 1. Axial profiles of uj^/uJc with defects: (a) defect location at sing^o = 1, (b) 
defect location at sing^o = — 1- 
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2.4- Even-parity eigenmode 

In contrast with figure 1(a), figure 1(b) illustrates a defect with tt < qz < 2tt that does 
not produce an odd-parity mode. Instead, an even-parity mode exists in this case. This 
mode is still described by (12) and (13), but the boundary condition at z = z now 
changes from (18) to 

~dE^(r;z) 



dz 



o, 



J Z=ZQ 



or, equivalent ly, 

F - Ge~ iqz ° = 0. 

The solvability conditions for (22), (23) and (28) now give 
cj — cjq e 



UJQ 



A 



-cos qz , 



sin qz , 



(27) 

(28) 

(29) 
(30) 



and we observe that A is positive (the wave field vanishes at infinity) when singz < 0, 
i. e. this even-parity eigenmode indeed requires tt < qz < 2tt for its existence. 



2. 5. Schrodinger equation for gap modes 

Equations (12) and (13) can be straightforwardly transformed into two independent 
second-order equations for F — G and F + G: 

16 d 1 d{F- G) 



dz ( u-uq _ e\ 
V 2 J 



dz 



+ {F-G) 



UJ 



16 d 



d{F + G) 



q 2 dz 



+ - 



dz 



+ (F + G) 



U 



CJ — CJQ 

CUq 



£ 
2 



0, 



(31) 



(32) 



Both of them can be further reduced to a time-independent Schrodinger equation when 
e is nearly constant. We assume e = (e) + u(z) with u(z) <C (e) and u(±oo) = 0. The 
discrete-spectrum modes in this case are very close to the tips of the gap. For the lower 
tip, we have 



CJQ 



2 



(33) 



with \S\ <C (e). We can then neglect u(z) and 5_ in the derivative term of (31) and get 



dHF-G) g*(e) 
dz^ + ~^ {F ~ G) 

Similar procedure for the upper tip gives 



S- + 



u(z) 



d 2 {F + G) , q 2 (e) 



dz 2 



+ 



16 



(F + G) 



-5 + + 



u(z) 



0. 



= 0. 



(34) 



(35) 



We note that any negative function u(z) acts as a potential well that supports spatially 
localised eigenmodes of (34) and (35) with 5_ > and S+ < 0, respectively. The 
eigenfrequencies of these modes belong to the spectral gap in the continuum. 
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3. Numerical implementation 

An ElectroMagnetic Solver (EMS) [T8] based on Maxwell's equations and a cold plasma 
dielectric tensor is employed to study the RLH spectral gap and gap eigenmode inside. 
The Maxwell's equations are expressed in the frequency domain: 



1 

Mo 



V xB = -zcjD+j a , 



(36) 
(37) 



where E and B are the wave electric and magnetic fields, respectively, D is the electric 
displacement vector, j a is the antenna current, and uj is the antenna driving frequency. 
These equations are Fourier transformed with respect to the azimuthal angle and then 
solved (for an azimuthal mode number m) by a finite difference scheme on a 2D domain 
(r; z). The computational domain is shown in figure 2. The quantities D and E are 



(r=0; z=0) 



Lz=20m 



10.047m 



VACUUM 



PLASMA 



0.213m 



Lr=0.2m 



J 

D; z = 0) > z 



Figure 2. Computational domain. The solid bar denotes a half-turn helical antenna. 
The dot-dashed line is the machine and coordinate system axis (r = 0). The coordinate 
system (r; (p; z) is right-handed with an azimuthal angle ip. 

linked via a dielectric tensor, [20j as follows: 

D = e [eE + ig(E x b) + (77 - e)(E ■ b)b], (38) 

where b = B o /5 is the unit vector along the static magnetic field and 

2 

<r — 1 _ 

UJ 2 

9 = ~ I] iff {u+ivly-ui a ' ( 39 ) 



*7 = l-£ 



a 



Oj(uO-\-ih>a) ' 



The subscript a labels particle species (electrons and ions); u) pa = ^n a q^/£om a is the 
plasma frequency, uj ca = q a B /m a is the cyclotron frequency, and v a a phenomenological 
collision frequency for each species. The static magnetic field is assumed to be 
axisymmetric with B 0r <C B 0z and B 0ip = 0. It is therefore appropriate to use a 
near axis expansion for the field, namely B 0z is only dependent on z and 

B 0ri r, ,) = ~\r™^>. (40) 
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A helical antenna is employed to excite an m = 1 mode in the plasma. The enclosing 
chamber is assumed to be ideally conducting so that the tangential components of E 
vanish at the chamber walls, i. e. 

E^{L r \ z) = E z (L r ; z) = 0, 

£ r (r;0)=^(r;0) = 0, (41) 

E r (r;L z ) = E (p (r ] L z ) = 0, 

with L r and L z the radius and the length of the chamber, respectively. Further, all field 
components must be regular on axis. 

4. Numerical results and discussion 

4-1. RLH mode in a straight cylinder 

We first use EMS to calculate the radial structure of the RLH mode in a straight cylinder 
with a uniform static magnetic field and recover the mode dispersion relation. We 
consider a single-ionised argon plasma with a radial density profile shown in figure 3(a). 
The surrounding vacuum chamber and the position of the rf-antenna are shown in 
figure 2. The static magnetic field in this simulation is B = 0.01 T and the electron- 
ion collision frequency on axis is h> ei \ r= $ = 4 x 10 6 s _1 . We run EMS for various 
frequencies of the antenna current, select the m — 1 azimuthal component of the 
plasma response and calculate the dominant axial wavenumber, k Z) of this component 
via Fourier decomposition. We observe an apparently preferred axial propagation 
direction for the m = 1 wave generated by the helical antenna. The collision frequency 
is sufficient to damp the radiated waves significantly towards the endplates but the 
dominant wavelength is still seen clearly in the plot of the excited field (figure 3(b)). 
The resulting relation between k z and uj is shown in figure 3(c). We note that k z scales 
as yJZj with frequency, in agreement with analytical scaling obtained in [19] . A least 
squares fit to the dispersion relation given by (17) with q/2 = k z shows that the form 
factor T equals 2.03 for the selected density profile. The radial profile of this m = 1 
mode is presented in figure 3(d). 

4.2. Spectral gap in a periodic structure 

In order to demonstrate the spectral gap, we run EMS for the same conditions as in 
section 14.11 except that the static magnetic field is weakly modulated along z and has 
the form B 0z = B (l + 0.25 cos qz) with B = 0.01 T and q = 40 m" 1 . The radial 
component of the static field is determined by (40). We scan the antenna frequency 
(for a fixed amplitude of the antenna current) and analyse wave propagation from the 
antenna. Figure 4 illustrates the results. Figure 4(a) is a plot of the on-axis signal at 
a 4.8 m distance from the antenna, and shows an interval of strong suppression (from 
36 Mrad/s to 44 Mrad/s), which represents the anticipated spectral gap. The location 
and width of this gap agree well with analytic estimates based on (16) and (17). Figure 
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Figure 3. Plasma density profile and RLH wave field in a straight cylinder: (a) radial 
profile of unperturbed plasma density, (b) azimuthal magnetic field of the m = 1 
mode for an illustrative frequency u = 31.4 Mrad/s (solid bar shows antenna location, 
and the endplates are located at z = m and z = 30 m respectively), (c) computed 
dispersion relation for the m = 1 mode (dots) and analytical scaling from (17) with 
q/2 = k z and T = 2.03 (dashed curve), (d) radial profile of the m = 1 mode at z = 20 m 
for uj = 31.4 Mrad/s. 

4(b) presents spatial distributions of wave energy for frequencies that are inside and 
outside of the spectral gap, from which we also see that wave propagation is evanescent 
inside the gap. 

4-3. Gap modes 

To form an eigenmode inside the spectral gap shown in figure 4, we introduce a defect in 
the otherwise periodic static magnetic field. Equations (25) and (29) suggest that this 
defect should be located at cos qz = 0, in order for the eigenmode frequency to be at 
the gap centre. The mode is then expected to have the shortest possible width (see (26) 
and (30)). Moreover, section [2731 and section 12741 indicate that the gap eigenmode can 
have either odd parity or even parity, depending on the defect profile. We will confirm 
this analysis by simulation results shown in figure 5 and figure 6. 

Figures 5(a)-(d) illustrate the odd-parity gap eigenmode associated with the defect 
shown in figure 1(a). In order to minimize the role of collisional dissipation and thereby 
obtain a sharper resonant peak in figure 5(c), we have reduced the collision frequency to 
0.1v e i in figure 5 simulation results. Figure 5(a) shows that the eigenmode is a standing 
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Figure 4. Identification of the spectral gap in simulations: (a) on-axis wave field at 
z = 15 m as a function of driving frequency; vertical lines show analytical predictions 
for the gap centre ujo and the upper and lower tips of the continuum cj + and c<j_, 
(b) spatial distributions of wave energy for antenna frequencies inside and outside the 
spectral gap (solid bar denotes the antenna location). 



wave localised around the defect. Its exponential envelope with a decay length of 1.25 m 
is consistent with analytical expectation from (26). The evanescent field of the mode 
is only weakly coupled to the distant antenna, but the resonance with the antenna 
frequency still allows the mode to be excited easily. The internal scale of this mode is 
close to twice the system's periodicity, consistent with Bragg's law. Figure 5(b) shows 
the on-axis profile of the mode electric field, implying an odd function of at zq. The 
corresponding magnetic field is, instead, an even function. Figure 5(d) presents a 
full view of the plasma response in [z\ uS) space with a clear eigenmode peak inside the 
spectral gap. 

We have also performed similar calculations for the defect shown in figure 1(b). 




Figure 5. Odd-parity gap eigenmode: (a) longitudinal profiles of the static magnetic 
field (solid line) and the rf magnetic field on axis for u = 40.2 Mrad/s (dots), together 
with the theoretically calculated envelope (dotted line), (b) longitudinal profile of 
(on- axis) at u = 40.2 Mrad/s (vertical dot-dashed line marks the defect location, and 
solid horizontal bar marks the antenna), (c) resonance in the dependence of the on-axis 
amplitude of the rf magnetic field on driving frequency at the location of the defect, 
(d) 3D plot of the on-axis wave field strength as a function of z and lj. 



This defect produces an even-parity mode seen in figures 6(a)-(d). Except for different 
symmetries, the features of the even and odd modes are apparently similar. 

The effects of collisionality on gap eigenmodes are shown in figure 7, from which 
we could see clearly a strength drop as the collisionality is increased. 

5. Summary 

We have shown that longitudinal modulation of the guiding magnetic field in a 
plasma column creates a spectral gap for radially localised helicon waves. Calculations 
performed with an EMS code reveal that this gap prohibits wave propagation along the 
column when the driving frequency of the rf antenna is in the forbidden range. The 
calculated width of the gap is consistent with analytical estimates. We have also shown 
that a discrete eigenmode can been formed inside the spectral gap by introducing a 
local defect to the periodic structure. Both the theoretical analysis and simulations 
demonstrate two types of gap eigenmode in the "imperfect" system: odd-parity mode 
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Figure 6. Even-parity gap eigenmode: (a) longitudinal profiles of the static magnetic 
field (solid line) and the rf magnetic field on axis for uj = 41.8 Mrad/s (dots), together 
with the theoretically calculated envelope (dotted line), (b) longitudinal profile of 
(on-axis) at uj — 41.8 Mrad/s (vertical dot-dashed line marks the defect location, and 
solid horizontal bar marks the antenna), (c) resonance in the dependence of the on-axis 
amplitude of the rf magnetic field on driving frequency near the location of the defect, 
(d) 3D plot of the on-axis wave field strength as a function of z and uj. 
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Figure 7. Effects of collisionality on gap eigenmodes: (a) results for odd-parity mode 
(figure 5(c)), (b) results for even-parity mode (figure 6(c)). 



and even-parity mode, depending the type of the defect employed. The gap mode is 
localised around the defect and represents a standing wave rather than travelling wave. 
Its distinctive feature is a resonant peak in the plasma response to the antenna current. 
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The gap eigenmode has two characteristic spatial scales: a short inner scale and a smooth 
envelope. The inner scale is nearly twice the system's periodicity, which is characteristic 
for Bragg's reflection; the envelope depends on the modulation amplitude and it scales 
roughly as the inverse width of the spectral gap. A plausible way to identify the gap 
mode in a linear device with multi-mirror configuration would be to use the end-plate of 
the machine as a controllable defect in the periodic system. This could make the mode 
observable with a modest number of mirrors in the machine. LAPD is an apparent 
candidate for such experiments, provided that dissipative processes in the plasma do 
not destroy the gap-mode resonance. 

Acknowledgments 

We would like to thank Dr. Guangye Chen for providing the EMS code and many 
instructions for its usage. This work was supported by the Chinese Scholarship 
Council [through scholarchip 2009611029], by the Postgraduate Research Award from 
the Australian Institute of Nuclear Science and Engineering, by the US Department 
of Energy under Contract No DE-FG02-04ER54742, and by the Australian Research 
Council through fellowship FT0991899. 

References 

[1] J.W. (Lord Rayleigh) Strutt. On the maintenance of vibrations by forces of double frequency, and 
on the propagation of waves through a medium endowed with periodic structure. Philosophical 
Magazine, 24:145-159, 1887. 

[2] P. W. Anderson. Absence of diffusion in certain random lattices. Phys. Rev., 109:1492-1505, Mar 
1958. 

[3] N. F. Mott. Conduction in non-crystalline systems .i. localized electronic states in disordered 

systems. PHILOSOPHICAL MAGAZINE, 17(150):1259, 1968. 
[4] Sajeev John. Strong localization of photons in certain disordered dielectric superlattices. Phys. 

Rev. Lett, 58:2486-2489, Jun 1987. 
[5] E. Yablonovitch, T. J. Gmitter, R. D. Meade, A. M. Rappe, K. D. Brommer, and J. D. 

Joannopoulos. Donor and acceptor modes in photonic band structure. Phys. Rev. Lett., 67:3380- 

3383, Dec 1991. 

[6] Alexander Figotin and Abel Klein. Localized classical waves created by defects. Journal of 

Statistical Physics, 86:165-177, 1997. 10.1007/BF02180202. 
[7] M. S. Chu, J. M. Greene, L. L. Lao, A. D. Turnbull, and M. S. Chance. A numerical study of 

the high-n shear Alfven spectrum gap and the high-n gap mode. Physics of Fluids B: Plasma 

Physics, 4(11):3713-3721, 1992. 
[8] D A D'Ippolito. Mode coupling in a toroidal sharp-boundary plasma, i. weak-coupling limit. 

Plasma Physics, 22(12):1091, 1980. 
[9] R. L. Dewar, R. C. Grimm, J. L. Johnson, E. A. Frieman, J. M. Greene, and P. H. Rutherford. 

Long- wavelength kink instabilities in low-pressure, uniform axial current, cylindrical plasmas 

with elliptic cross sections. Physics of Fluids, 17(5):930-938, 1974. 
[10] R. Betti and J. P. Freidberg. Stability of Alfven gap modes in burning plasmas. Physics of Fluids 

B: Plasma Physics, 4(6):1465-1474, 1992. 
[11] N. Nakajima, C. Z. Cheng, and M. Okamoto. High-n helicity-induced shear Alfven eigenmodes. 

Physics of Fluids B: Plasma Physics, 4(5):1115-1121, 1992. 



Gap eigenmode of radially localised helicon waves 



14 



[12] Ya. I. Kolesnichenko, V. V. Lutsenko, H. Wobig, Yu. V. Yakovenko, and O. P. Fesenyuk. Alfven 
continuum and high-frequency eigenmodes in optimized stellar ators. Physics of Plasmas, 
8(2):491-509, 2001. 

[13] H.H. Duong, W.W. Heidbrink, E.J. Strait, T.W. Petrie, R. Lee, R.A. Moyer, and J.G. Watkins. 

Loss of energetic beam ions during tae instabilities. Nuclear Fusion, 33(5) :749, 1993. 
[14] R. B. White, E. Fredrickson, D. Darrow, M. Zarnstorff, R. Wilson, S. Zweben, K. Hill, Yang Chen, 

and Guoyong Fu. Toroidal Alfven eigenmode-induced ripple trapping. Physics of Plasmas, 

2(8):2871-2873, 1995. 

[15] C.Z. Cheng L. Chen and M. S. Chance. High-n ideal and resistive shear Alfven waves in tokamaks. 

Annals of Physics, 161(1):21, 1985. 
[16] King-Lap Wong. A review of Alfven eigenmode observations in toroidal plasmas. Plasma Physics 

and Controlled Fusion, 41(1):R1, 1999. 
[17] Yang Zhang, W. W. Heidbrink, H. Boehmer, R. McWilliams, Guangye Chen, B. N. Breizman, 

S. Vincena, T. Carter, D. Leneman, W. Gekelman, P. Pribyl, and B. Brugman. Spectral gap of 

shear Alfven waves in a periodic array of magnetic mirrors. Physics of Plasmas, 15(1):012103, 

2008. 

[18] Guangye Chen, Alexey V. Arefiev, Roger D. Bengtson, Boris N. Breizman, Charles A. Lee, and 

Laxminarayan L. Raja. Resonant power absorption in helicon plasma sources. Physics of 

Plasmas, 13(12):123507, 2006. 
[19] Boris N. Breizman and Alexey V. Arefiev. Radially localized helicon modes in nonuniform plasma. 

Phys. Rev. Lett, 84:3863-3866, Apr 2000. 
[20] V. L. Ginzburg. The Propagation of Electromagnetic Waves in Plasmas. Pergamon, Oxford, p. 

84, 1964. 



